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We study the main spectral features of the gamma-ray fluxes observed by the High Energy Stereo- 
scopic System (HESS) from the J1745-290 Galactic Center source during the years 2004, 2005 and 
2006. In particular, we show that these data are well fitted as the secondary gamma-rays photons 
generated from dark matter annihilating into Standard Model particles in combination with a sim- 
ple power law background. We present explicit analyses for annihilation in a single standard model 
particle-antiparticle pair. In this case, the best fits are obtained for the uu and dd quark channels 
and for the W + W~ and ZZ gauge bosons, with background spectral index compatible with the 
Fermi-Large Area Telescope (LAT) data from the same region. The fits return a heavy WIMP, with 
a mass above ~ 10 TeV, but well below the unitarity limit for thermal relic annihilation. 

PACS numbers: 04.50.Kd, 95.36.+X, 98.80.-k 

I. INTRODUCTION 

Numerous and different data collected during the last years have established the Standard Cosmological Model as 
a simple framework showing a remarkable agreement with observations. This model is based on Einstein's General 
Relativity and a homogeneous and isotropic ansatz for the metric. However, in addition to ordinary matter, two 
new elements need to be added: a cosmological constant or other form of dark energy to account for the late time 
acceleration of the universe, and dark matter (DM) to explain the formation and dynamics of cosmic structures. 
However, despite the multiple efforts, the fundamental nature of DM remains still as an open problem. There are 
strong astrophysical evidences for DM from galactic to cosmological scales, but the interactions with ordinary matter 
have not been probed beyond gravitational effects. In this sense, both direct and indirect DM searches are fundamental 
to explore particle models of DM. 

In the framework of indirect searches, the observation of gamma-ray fluxes from astrophysical sources represents 
one of the main approaches to the DM puzzle. If DM annihilates into SM particles, the subsequent chains of decay and 
hadronization of unstable products produce gamma-ray photons generically. The observation of this signal is highly 
affected by astrophysical uncertainties in the gamma-ray backgrounds and in the DM densities and distribution. 
On the other hand, depending on both astrophysical and particle physics models of DM, the signature could be 
distinguishable from the background. Appealing targets for gamma-ray observations of annihilating DM arc mainly 
selected according to the abundance of DM in the source and their distance. Galaxy clusters, dwarfs spheroidal 
galaxies (dSph) or the galactic center (GC) of the Milky Way are traditional targets of interest. Galaxy clusters 
are very rich in DM, but they are very distant objects. DSphs are also characterized by high DM densities but at 
much shorter distances. In any case, despite their proximity, observations of dSphs have not been able to detect any 
gamma-ray flux signature over the background so far [TJH]- On the other hand, the GC represents a very close target, 
but the complexity of the region, due to the high concentration of sources, makes the analysis quite involved. 

In this work, we will focus on the analysis of very high energy (VHE) gamma-rays coming from the GC. Different 
observations from the GC have been reported by several collaborations such as CANGAROO [3 , VERITAS [4 , HESS 
[SJE], MAGIC [7] and Fermi-LAT [8, 9 . In particular, we will analyze the data collected by the HESS collaboration 
during the years 2004, 2005, and 2006 associated with the HESS J 1745-290 GC source [B]. The interpretation of 
these fluxes as DM signal has been widely discussed in the literature from the very early days of the publication of 
the observed data [T0"l - fT3] . It was concluded that the spectral features of these gamma-rays disfavored the DM origin 
[T^llbl] . However, in recent studies |151I16| . it has been pointed out that these fluxes are compatible with the spectral 
signal of DM annihilation provided it is complemented with a simple background. This extra source of gamma-rays 
can be originated by radiative processes generated by particle acceleration in the vicinity of Sgr A East supernova 
and the supermassive black hole [17J. 

In this work, we present a systematic study of this assumption: In Section ^ we show an analysis of the source, 
while a brief review of the gamma-ray flux coming from annihilating DM in galactic sources is presented in Section 
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|III| The fit of the HESS data with a total fitting function given by the combination of the background power law 



component with annihilating DM signature is discussed in Section IV In Section V, we include in the analysis the data 
collected by Fermi-LAT from the same region. We summarize the main results of our work and prospects for future 
analyses in Section |VI| Finally, an appendix provides useful details for reproducing the statistical study performed in 
these analyses. 

II. HESS J1745-290 DATA 

DM is expected to be clumped in the central region of standard galaxies. In this sense, the central part of the Milky 
Way could be an important source of gamma-rays produced in the DM annihilation processes. Because of its closeness, 
the GC represents a very appealing target for the indirect search of DM, but the complex nature of this area makes 
the identification of the sources quite difficult. Several sources have been detected not only in the gamma-ray, but also 
in the infrared and X-ray ranges of the spectrum. The absence of variability in the collection of HESS J 1745-290 data 
in the TeV scale, during the years 2004, 2005 and 2006, suggests that the emission mechanism and emission regions 
differ from those invoked in the variable infrared and X-ray emissions |18j . The significance of the signal reduces to 
few tenths of degree [5] , but the angular distribution of the very high energy gamma-ray emission shows the presence 
of an adjunctive diffuse component. The fundamental nature of this source is still unclear. These gamma-rays could 
have been originated by particle propagation j5J [T7] in the neighborhood of the Sgr A East supernova remnant and 
the supermassive black hole Sgr A, both located at the central region of our galaxy [l"4"l |19j . If it was originated by 
the DM distribution, the morphology of the source requires a very compressed DM structure. In fact, it has been 
claimed [20]. [21] that baryonic gas could account for similar effects when falling to the central region, modifying the 
gravitational potential and increasing the DM density in the center (see however [22, 23J). If this is the case, the DM 
annihilating fluxes are expected to be enhanced by an important factor with respect to DM alone simulations, such 
as the classical NFW profile [21] . 

In previous works [6] [19], important deviations from a power law spectrum have been already reported, and a cut-off 
at several tens of TeVs has been proved as a distinctive feature of the spectrum. For instance, the observed data were 
compared to a simple power law: 

^ = *o-(— V 1 " (1) 
dE \GeVj ' 1 ' 

and a power law with high energy exponential cut-off: 



-r 



dE -MGevJ (2) 

where $ is the flux normalization, Y is the spectral index and E cut is the cut-off energy. The measured spectrum for 
the three-year dataset ranges from 260 GeV to 60 TeV. We have reproduced these analyses with the results that can 
be found in Figs. [I] and [2] They are consistent with previous works [5] HH] , since the spectrum is well described by a 
power law with exponential cut-off (Fig. |2j, while a simple power law is clearly inconsistent (Fig. [TJ. 

III. GAMMA-RAYS FROM DARK MATTER ANNIHILATION 

If the signal observed from the GC is a combination of gamma-ray photons from annihilating DM and a background, 
the total fitting function for the observed differential gamma-ray flux will be: 

d$Tot _ d$> B g &DM 

dE ~ dE dE ■ U 
We will assume a simple power law background parameterized as: 

On the other hand, the differential gamma-ray flux originated from DM annihilation in galactic structures and 
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Figure 1: Collection of the HESS data (2004,2005 and 2006) fitted with two simple power law background. The full line (Power law 1) 
shows the fit with 23 dof, without take into account the upper limit constraints on the flux given by the last three points. They are 

included in the Power law 2 fit with 26 dof (dotted line). The parameters of the fit can be found in the table. Because of the poorness of 
both the fits (x 2 /dof = 2.48 with V = 2.32 ± 0.03 for Power law f , and x 2 /dof = 5.30 with T = 2.43 ± 0.03 for Power law 2), they 
represent just an upper limit of the flux. These results are consistent with previous studies l6l 1191 . 
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Figure 2: The same as in Fig. ^ but the power law is modified with an exponential energy cut-off. The signal is well fitted in this 

case with the results given in the attached table. 



substructures can be written generically as: 

dE 



channels 

E 



dE 



A^(J) AQ 
4nM 
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where (o^d) are the thermal averaged annihilation cross-sections of two DM particles into SM particles (also labeled 
by the subindex i). We are assuming that the DM particle is its own antiparticle. M is the mass of the DM particle, 
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and the number of photons produced in each annihilating channel dN-JdE. involves decays and/or hadronization of 
unstable products such as quarks and leptons. Because of the non-perturbative QCD effects, the calculation of dNi/dE 
requires Monte Carlo events generators such as PYTHIA [24 . In our analysis, we will focus on gamma-rays coming 
from external bremsstrahlung and fragmentation of SM particle-antiparticle pairs produced by DM annihilation. We 
will ignore DM decays, the possible production of monoenergetic photons, N-body annihilations (with N > 2), or 
photons produced from internal bremsstrahlung, that are model dependent. In particular and in order to simplify 
the discussion and provide useful information for a general analysis, we will consider DM annihilation into each single 
channel of SM particle-antiparticle pairs, i.e. 



~dE~ ~ Ai '^E' (6) 

where 

2 _ (aiv)AQ{J) AQ 

is a new constant to be fitted together with the DM particle mass M, and the background parameters B and T. The 
astrophysical factor (J) An of Eq. ^ will be also indirectly fitted by means of the parameter Ai. In its most general 
expression, it is computed in the direction W defined by the line of observation towards the GC: 



1 



An 



AO 



(J)An = l^ dfi p 2 [r(l)]dl(V) , (8) 



where I is the distance from the Sun to any point in the halo. The radial distance r is measured from the GC, and 



is related to I by r 2 — I 2 + Dq — 2DqIcos^, where Dq ~ 8.5 kpc is the distance from the Sun to the center of the 



Galaxy. The distance from the Sun to the edge of the halo in the direction 6 is l max = Dq cos 9 + yr 2 — D & sin 0. 

The photon flux must be averaged over the solid angle of the detector, that is typically of order Aft — 2tt(1~ cos ~ 
10~ 5 for detectors with sensitivities in the TeV energy scale, such as the HESS Cherenkov telescopes array. The dark 
halo in the GC is usually modeled by the NFW profile [25 : 

TV 

r[r — r s y 

where N is the overall normalization and r s the scale radius. This profile is in good agreement with non-baryonic 
cold DM simulations of the GC. In this case and accounting for just annihilating DM, the astrophysical factor is: 
(J(T) W ) — • 10 23 GeV 2 cm -5 , value that we will use as standard reference. 



IV. SINGLE-CHANNEL FITS 



As commented before, the particle model part of the differential gamma-ray flux expected from the GC is simulated 
by means of Monte Carlo event generators, such as PYTHIA [24J. However, the fact that simulations have to be 
performed for fixed DM mass implies that we cannot obtain explicit M dependence for the photon spectra. In order 
to overcome this limitation, the simulated spectra of each annihilation channel has been reproduced with the analytic 
fitting functions provided in Ref. |26j in terms of the WIMP mass, by means of mass dependent parameters. The 
combination of such simulated spectra with a power law background (Eq. (J3j) ) is finally fitted. We assume a typical 
experimental resolution of 15% (AE/E ~ 0.15) and a perfect detector efficiency. 

For quarks (except the top) electron and r leptons, the most general formula needed to reproduce the behaviour of 
the differential number of photons in an energy range may be written as: 
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The value of the different parameters change with the SM particle annihilation channel and in some cases, with 
the range of the WIMP mass. The cases of interest are described below and the value of the parameters reported in 
Appendix A. In the case of the electron-positron channel, the only contribution to the gamma-rays flux comes from the 
bremsstrahlung of the final particles. Therefore, in the previous expression ( 10 1 the exponential contribution is absent, 

<1 = "qed/tt, and p = (M/m e 
approximation (Fig. [3|. 



) . This choice of the parameters corresponds to the well-known Weizsacker- Williams 



In the case of the /i^/i channel, the exponential contribution in the expression above ( 10 ) is also absent. A proper 
fitting function for such a channel can be written as: 
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All the parameters are here mass dependent and their expression for a range of mass 10 3 GeV < M < 5 x 10 4 
are reported in Tab. |T] The best fit for the a lepton is shown in Fig. [4] 



(11) 



GeV 



The t + t~ spectrum needs the entire Eq. ( 10 1 for an accurate analysis. The value of the mass independent parameter 
and the expression of the mass dependent ones used in this work are reported in Tab. [IT] The best fit for the r lepton 
is shown in Fig. [5j 

The value of the mass independent parameter and the expression of the mass dependent ones in ( 10 1 for the uu 
channel are reported in Tab. 



III| The analytic fitting function is given for mass values up 8000 GeV, because of 
limitations in the Monte Carlo event generator PYTHIA software. An extrapolation up to larger values of the mass 
has been performed. In this case the best fit is shown in Fig. [6] The parameters for the dd channel are given in Tab. 
IV This channel provides the best fit (See Fig. [7]) of all those considered in the paper and is used as a reference for 



comparison with other channels. The parameters for the ss, cc and bb are reported in Tabs, 
fit of these hadronic channels are shown in Figs. [8] [9] and 10 
The ti needs a different fitting function: 
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The value of the mass independent parameter and the expression of the mass dependent ones (in the selected range 



of mass) are reported in Tab. VIII The best fit for the tt channel by using Eq. (12 1 is shown in Fig. 11 
For the W and Z gauge bosons the parametrization is: 
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IX The best fits for the W+W~ and ZZ 



where the value of the parameters used in this study are reported in Tab 
channels are shown in Figs. [l2]and[T3j 

The best fit is provided by the dd channel with x 2 /dof = 0.73 for a total of 24 degrees of freedom (dof). In any 
case, other hadronic channels such as uu (see Fi g. [6| ) or ss, also provide very good fits within la. In the same way, 

W+W 



121 or ZZ (Fig. 



13| channels are consistent with 
or but also 



softer spectra as the one provided by tt (see Fig. IllE W + W~ (Fig. 

data without statistical significance difference. On the contrary, leptonic channels (not only e n 

t + t~, Figs. [3] [4] and [5]), cc and bb channels (Fig. [9] and 10) are ruled out with more than 99% confidence level when 
compared to the best channel. It is interesting to note that taking into account all the channels that provide a good 
fit, the DM mass is constrained to 15 TeV < M < 110 TeV within 2a. The lighter values are consistent with hadronic 
annihilations (uu) and the heavier ones with the annihilation in ti, that is more similar to electroweak channels. 



FERMI 1FGL J1745. 6-2900 DATA 



It has been argued that the Fermi-LAT source 1FGL J1745. 6-2900 and the previously considered HESS source 
J 1745-290 are spatially coincident j^Tj. In |H| data from the first 25 months of observations of the mentioned Fermi- 
LAT source have been analyzed. It has been shown that the observed spectrum in the range 100 MeV- 300 GeV can 
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Figure 3: M(TeV), A(10- 7 cm-'s" 1 / 2 ), B(10" 4 GeV" 1 / 2 ^-^). Best fit to the HESS J1745-290 collection of data (years 2004,2005, 
and 2006 [6]) in the case that the DM contribution came entirely from annihilation into e + e — . The full line shows the total fitting 
function. The poor quality of the fit is evident (x 2 /dof = 2.09). The dotted line shows the simple power law background component 
with spectral index T = 2.55 ± 0.06, while the DM component with a 15% of resolution uncertainity (R. U.) is given by the double-dotted 
line with a normalization parameter A. The dotted-dashed line shows the contribution of annihilating DM into electron-positron pairs 
without taking into account such a resolution. On the top right table, we present a summary of the best fitting parameters. A% 2 
correponds to the difference with respect to the best fit channel (dd). On the bottom right table, we show the correlation matrix 

elements of this statistical analysis. 
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Figure 4: The same as Fig. [3] in the case of the /x + /i annihilation channel. The fit is only slightly better than for the 



be very well described by a broken power law with a break energy of E^ r — 2.0_"'q GeV, and slopes 1\ = 2.20 ± 0.04 
and T2 = 2.68 ± 0.05 for lower and higher energies than E\> r respectively. Notice that the fitted value of from 
Fermi-LAT data is in very good agreement with the spectral index of the diffuse background obtained from HESS 
data in our previous analyses. Indeed, at the 95% confidence level, the allowed range for the spectral index of the 
diffuse background is 2.6 < T < 3.7. In this case, the lower values are consistent with all the allowed channels, but 
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Figure 5: The same as Fig. |3]in the case of the t+t - annihilation channel. This channel provides the best leptonic fit, although it is 

not enough to be consistent with data. 




the higher values are only accessible to the light quark channels. In any case all the channels that provide a good fit 
to HESS data are also consistent with the spectral index observed by Fermi-LAT. 



In Fig. 14 we show the case of the W + W~ channel to illustrate this consistency. Both the signal and background 
parameters are compatible with the W + W~ channel fit without the Fermi-LAT data. With the new data, the analysis 
even improves to x 2 /dof — 0.75. This interpretation implies that the Fermi-LAT telescope is able to detect just the 
background component of the total energy spectrum of the gamma-ray emission associated with DM. 

With the fit of the parameter A (Eq. Q) and by assuming a standard thermal cross section of (av) = 3 • 
10~ 26 cm 3 s -1 , we can get an estimation of the astrophysical factor (J). We find 10 25 GeV 2 cm~ 5 < (J NFW ) < 
10 26 GeV 2 cm~ 5 . It implies that the boost factors b = (J)/(J NFW ) spread on a range between two and three orders 
of magnitude. It is interesting to note that the enhancement of the distribution of DM required to fit the data is 
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Figure 7: The same as Fig. [3] in the case of the dd annihilation channel. The x 2 /dof value associated with this analysis is the lowest 
one obtained in this work and it has been taken as reference to compute Ax 2 . 
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Figure 8: The same as Fig. [SJin the case of the ss annihilation channel. 

compatible with the expectation of N-body simulations when the effect of the baryonic matter in the inner part of 
GC are taken into account |21j . 

VI. CONCLUSIONS 

In this work, we have analyzed the possibility of explaining the gamma-ray data observed by HESS from the central 
part of our galaxy by being partially produced by DM annihilation. The complexity of the region and the ambiguous 
localisation of the numerous emitting sources inside this region, seems to validate the hypothesis of a background 
component. We have proved that even DM annihilations into single channels of the SM provide good fits provided that 
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Figure 9: The same as Fig. [3] but for the cc channel. This is the quark channel most disfavored by the HESS data. 
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Figure 10: The same as Fig. [3] for the I 
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b annihilation channel. Together with the cc channel, is the only hadronic channel disfavored 
by the HESS data.) 



the DM signal is complemented with a diffuse background compatible with Fermi LAT observations. The fits returns 
a DM mass between 15 TeV < M < 110 TeV and a background spectral index between 2.6 < V < 3.7 within 2a. 
Some channels are clearly preferred with respect to others, such as the dd for the quark channels with x 2 1 dof = 0.73, 
or the W + W~ and ZZ channels with x 2 /dof = 0.84 and x 2 /dof = 0.85. On the contrary, the leptonic channels are 
seriously disfavored. The morphology of the signal is consistent with compressed dark halos by taking into account 
baryonic dissipation |20[ 121] with boost factors of b ~ 10 3 for typical thermal cross sections. 

The DM particle that may have originated these data needs to be heavier than ~ 10 TeV. These large DM masses 
required for fitting the HESS data are not in contradiction with the unitarity limit for thermal relic annihilation. For 
example, this limit is around 110 TeV for scalar DM particles annihilating in s-wave (for JIdm^ 2 = 0.11) [151 12"8"] . 

On the contrary, these heavy DM particles are practically unconstrained by direct detection experiments or particle 



10 



CM 10' 



O 

> 



CD 

X -1 
3 10 1 



LU 



10"' 




Channel 



tt 



HESS data (2004,2005,2006) 
Background 
Total fit 

t channel with 15% R.U. 
t channel with 0% R.U. 



M 

A 
B 

r 

X 2 /dof 

Ax 2 
6 



87.7 ±8.2 
3.68 ±0.34 
6.07 ±3.34 
2.86 ±0.19 

0.88 

3.6 
3116 ±820 





tt 


M 


A 


b r 


M 


1 






A 


-0.748 


1 




D 


-0.175 


0.721 


i 


r 


-0.216 


0.756 


0.998 1 



E (GeV) 



Figure 11: The same as Fig. [3] in the case of the tt annihilation channel. The qualitative features of this channel are more similar to 
the electroweak annihilations (see following figures) than to the hadronic ones (as it is shown in the previous figures). The quality of the 

fit is also good. 
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Figure 12: The same as Fig. plfor the W + W annihilation channel. 



colliders An interesting example of such type of DM candidates which could have high enough mass and account 
for the right amount of DM in the form of a thermal relic, is the branon |30U31j . Branons are new degrees of freedom 
corresponding to brane fluctuations in brane-world models. In general, they are natural candidates for DM because 
they are massive fields weakly interacting with the SM particles, since their interaction is suppressed with the fourth 
power of the brane tension scale /. For masses over 100 GeV, the main contribution to the photon spectra comes 
from branons annihilating into gauge bosons ZZ, and W + W~ . In [TS] it was shown that a branon DM with mass 
M ~ 50.6 TeV, provides an excellent fit to HESS data. The corresponding background being compatible with Fermi- 
LAT data. The compatibility of its thermal abundance with the WMAP constraints |32_ , demands a cross section of 
(<jv) — (1.14 ± 0.19) ■ 10~ 26 cm 3 s~ 1 , what is equivalent to a brane tension of / ~ 27.5 TeV. 
The analysis of other cosmic rays |33j from the GC and from other astrophysical objects is fundamental to cross 
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Figure 13: The same as Fig. [3]in the case of the ZZ annihilation. The results for the ZZ channel are very similar to the W + W one. 

Both electroweak channels are consistent with the data. 




Figure 14: The same as in Fig. [l2]but taking into account the Fermi-LAT data between 2 GeV < E < 100 GeV [9], that extends the 
dof to 27. The background parameters fitted by these data are A = (4.44 ± 0.34) ■ 10~ 'cm - 1 s~ 1 / 2 and V = 2.63 ± 0.02, in agreement 
with the spectral index for this spectral region already found when the Fermi-LAT data are fitted with a broken power law [5]. The 
quality of the fit is evident with a x 2 /dof = 0.75. From this interpretation of the data, Fermi-LAT instruments are sensitive just to the 

background component of the signal. 



check the hypotheses considered in this work. Updated analyses of this kind of signals for heavy DM combined with 
simple background components would be of great interest in order to constrain the possible DM origin of the studied 
gamma-ray fluxes. 
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APPENDIX A: Fitting function parameters 

In the following tables, we show explicitly the parameters used in the fitting functions (read [26 for further details). 



Parameter 


Fitting power law(s) 


P 


4530Af ubba 


q 


0.00230Af" uy11 + 0.00291Af ou;:i48 


I 


0.626M- UU3UU + 16.4M" 1 - 34 



Table I: Mass dependent parameters for the fitting function of the photon spectrum coming from the fj, + /i channel. 



Parameter 


Fitting power law 


ni 


-7.00M" iy9 + 119M- (,lbA + 9.09 


P 


3.07M 155 



Channel 


fll 


h 


b 2 


n 2 




di 


C2 


d 2 


q 


T 


14.7 


5.40 


5.31 


1.40 


2.54 


0.295 


0.373 


0.470 


0.00260 



Table II: Mass dependent and independent parameters for the r + r channel fitting function. 



Parameter 


Fitting power law 


bi 


2.Q6M 1 ™ 


ni 


2.91M- uabi + 1.90Af uuiY ^ 


ri2 


0.0587M°- 14b ' + 0.848M" U - 145 


di 


0.317M- uuauu + 0.403M- uabi 


P 


4 74M U.839 



Channel 


Oi 


b 2 


Cl 


c 2 


q 


uu 


5.58 


5.50 


0.315 


0.0 


9.30 • 10~ 4 



Table III: Parameters for the uu channel fitting function. 



Parameter 



Fitting power law 



0.0485 



3.39Af 



-0.!)!)3 



-0.00113 



"1 



21.8M 



+ 2.25M 



-0219 



0.0573 



"2 



0.848Af 



0.161M 



ci 0.722M-° 27U + 0.0544A/ 



OTTTT 



JT72W 



0.168M 



Channel 


a\ 


b 2 


di 


C2 


q 


dd 


5.20 


5.10 


0.410 


0.0260 


1.40 • 10" 4 



Table IV: The analogous set of parameters shown in Tab. Ill but for the dd channel 
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Parameter 


Fitting power law 


bi 


4. 54M (™ 


n 2 


3.68M" iUi + 0.744M- uuaM 


di 


0.621M" u - bY4 + 0.414M" UU!38S 


V 


12.8M"- 732 



Channel 




m 


b 2 


Cl 


9 


ss 


4.83 


2.03 


6.50 


0.335 


g = 2.40 ■ 10~ 4 



Table V: Parameters for the fitting function of the photon spectrum coming from DM annihilation in ss channel. 



Parameter 


Fitting power law 


h 


9.90M- Uiau 


ni 


4.14M- u - i4 » 


Cl 


0.210M 00951 


di 


1.50Af- Ub3 ^ + 0.479M- uuy4 ^ 


P 


8.11M"- 812 



Channel a\ b 2 n 2 c 2 



5.58 7.90 0.686 0.0 9.00 • 10~ 4 



Table VI: The same as in Tab. Ill for the cc channel 



Parameter 


Fitting power law 


bi 


152M- u - 4b2 


m 


18.7M- UMti 


n 2 


0.707M" UUi;iy 


Cl 


0.328M UU44Y 


di 


0.449M" UU5M 


P 


11.8M ue4i 



Channel 


at 


b 2 


C2 


d 2 


q 


bb 


10.0 


11.0 


0.0151 


0.550 


2.60- 10" 4 



Table VII: Parameters for the bb channel fitting function. 



Parameter 


Fitting power law(s) 


h 


16.4M- U - U4UU 


m 


0.559M" UU3Ya 


c 2 


8910A/- 3 - 23 


p 


5.78 • 10- 5 M 189 


q 


0.133M U488 


i 





Channel 


a% 


Cl 


di 


d 2 


tt 


290 


1.61 


0.19 


0.845 



Table VIII: Mass depend and independent parameter of the fitting function given by Eq. ( 12 I for the tt channel. 
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Channel 


ai 


ni 


9 


W 


25.8 


0.510 


3.00 


z 


25.8 


0.5 


3.87 



Channel 


h 


Cl 


di 








9 .29M- UU139 


0.743M UU33i 


0.265M- UUi;;iY 


10 5 M~ ii3 + 285M U - U794 


0.943M°- ,)u8b2 


ZZ 


936M -u.uoviu 


0765M U.UU98U 


0.272M- UUU99U 


85505M-°- ibb + 0.476M U9b4 


0.884M"- UiYt> 



Table IX: Parameters for the fitting functions of the electroweak gauge bosons channels. 



